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Abstract 

We compare several definitions for number-conserving cellular au- 
tomata that we prove to be equivalent. A necessary and sufficient 
condition for cellular automata to be number-conserving is proved. 
Using this condition, we give a linear-time algorithm to decide number- 
conservation. 

The dynamical behavior of number-conserving cellular automata is 
studied and a classification that focuses on chaoticity is given. 

Keywords: Cellular automata, decidability, discrete dynamical systems 



1 Introduction 

Cellular automata are widely known as formal models for complex systems 
ruled by local interactions [28, 25, [H| 23, ^6|. Despite of their apparent 



simplicity they display a richness of behaviors most of which are not yet fully 
understood. The classification of these behaviors has a main obstacle: most 
of the interesting properties about them are undecidable |flq| . Moreover, 
some other properties are "dimension sensitive"; for example, surjectivity is 
decidable in polynomial time for one-dimensional CA and undecidable for 



dimension > 2 [15, 10 



We prove that number-conservation, a widely used property, can be de- 
cided in linear time (in the size of the CA) in any dimension. More precisely, 
it has linear time complexity if we assume that the access to the transition 
table costs constant time, independently of the dimension of the space. Oth- 
erwise the complexity is ©(slogs), where s is the size of the transition table. 
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A preliminary work had been carried out for one- dimensional CA on circu- 
lar configurations (periodic boundary conditions) by Boccara et al. in @, ||]. 
The decision problem in higher dimensions and general boundary conditions 
that we solved in the present paper is rather surprising and could not be 
expected from earlier works. 

CA that satisfy number-conservation property are called number- conser- 
ving (NCA); roughly speaking they conserve the number of bits during their 
evolution. This class of CA has been introduced by [24] and is receiving a 
growing interest in Physics as a source for models of particles systems ruled 
by conservations laws of mass or energy. 

Great attention has received the application of NCA to practical situa- 
tions such as highway traffic flows. In this context different definitions of 
"number-conservation" were given according to the type of boundary con- 
ditions assumed (periodic, finite, infinite) |], ||. In Section || we prove 
that all of them are equivalent. We also give another equivalent definition 
in terms of densities. Section ||| is devoted to the proof of our main re- 
sult: the linear time decision algorithm. In the second part of the paper we 
study NCA dynamics with focus on chaotic behavior. Number-conservation 
imposes strong "simplifications" in the dynamics allowing (in most cases) 
clear-cut situations that are not observed in the general case. 

We stress that the class of NCA is not only big enough to give reliable 
models for car traffic and similar models allowing the study of traffic jams, 
critical traffic masses, crossing-roads, car accidents etc. But also they are 
interesting from a pure computer science point of view. For instance Morita 
et al. have given a reversible two-dimensional NCA capable of universal com- 
putation pOl |2l. 



2 Definitions 

We present in this section several definitions of number- conserving cellu- 
lar automata that we prove to be equivalent. Proofs are given in a one- 
dimensional space but can be generalized to any dimension. 

Cellular automata are formally defined as quadruples (d,S, N,f). The 
integer d is the dimension of the space the cellular automaton will work 
on. S = {0,1,... , s} is called the set of states. The neighborhood N = 
(xi,...x v ) is a v-tuple of distinct vectors of Z rf . The ccj's are the relative 
positions of the neighbor cells with respect to the cell, the new state of which 
is being computed. The states of these neighbors are used to compute the 
new state of the center cell. The local function of the cellular automaton 
/ : S v i— > S gives the local transition rule. 

A configuration is an application from 7L d to S. The set of configurations 
is S z<i on which the global function A of the cellular automaton is defined 



2 



via /: 



Vc G 5 zd ,Vi G Z rf ,A(c)(i) = /(c(i + a;i),... ,c(i + x„)). 

We denote by Cp the set of (spatial) periodic configurations, (i.e., con- 
figurations that are periodic in all dimensions of the space). We denote by 
7r(c) the period of a periodic configuration (i.e., a vector.) The expression 
< k < 7r(c), k € Z d means VI < i < d, < k{ < 7r(c)j where d is the 
dimension of the space and ki is the i-th component of k. 

Definition 1 Let A be a d- dimensional cellular automaton. A is said to be 
periodic-number-conserving (PNC) iff 

0<fc<7r(c) 0<fc<7r(c) 

We denote by Cp the set of finite configurations, (i.e., those configurations 
that are almost everywhere equal to zero — non-zero in a finite number of 
cells). 

Definition 2 Let A be a d-dimensional cellular automaton. A is said to be 
finite-number-conserving (FNC) iff 



VcGC F , ^2c(i) = ^2A(c)(i 



We denote by C the set of all configurations. 

Definition 3 Let A be a d-dimensional cellular automaton. A window is 
a hypercube of Z d centered in thus determined by its size. Consider the 
sequence of windows {F n } of size 2n + 1 and denote by fj, n (c) the sum of 
states in F n of a configuration c G C. A is said to be number- conserving 
(NC) iff 

Vc G C, lim -7-77-^rT = 1. 

n^oo n n (A(c)) 

Remark that all these definitions imply that /(0, 0, . . . , 0) = 0, i.e., is 
a quiescent state. The configuration with all quiescent states is denoted uz. 

We are going to prove that all previous definitions are equivalent. The 
proof is given in dimension 1; the generalization to higher dimensions is 
straightforward . 

Proposition 1 A cellular automaton A is FNC if and only if it is PNC. 
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Proof. We first prove that if A is PNC then it is FNC. We show that if A is 
not FNC then it is not PNC either. 
Let A be not FNC, then 

3c £ Cp, 3a, /3 G N, a = ^ c(i) / ^ A(c)(i) = /?. 

Let be the length of the non-zero part of c and Z the radius of the 
neighborhood of A. Then, the length of the non-zero part of A(c) is at most 
k + 21 (see Figure ||) . 

Let c be the periodic configuration constructed from c in the following 
way. We add I zeros at both sides of the non-zero part of c and we construct 
a periodic configuration with the segment obtained (such a configuration has 
a period p = k + 21, see Figure |l|). 
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Figure 1: Construction of the periodic configuration 
We have 

0<i<p 0<i<p 

hence A is not PNC. 

We now prove that if A is FNC then it is PNC. We show that if A is not 
PNC then it is not FNC either. 

Let I be the radius of the neighborhood. As A is not PNC, 

3p,a,f3 E N,3c G C P , ^ c(i) = a ^ (3 = ^ A (c)(i), 

0<i<p 0<i<p 

where p is the period of c. 
We denote by 

c p the p-length segment of c centered at the origin; 
CpSuffix; the l-length suffix of c p ; 
CpPrefixi the Z-length prefix of c p . 
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We construct a finite configuration c G Cf a as shown in Figure ||: we 
repeat c p x times and we complete the segment obtained in this way with 
CpPrefixj at the right-hand-side and with CpSuffix; at the left-hand-side of the 
segment. 



A(c): I 



pSuffix, 



pPrefix. 







x times 
A(h, A(c) p 



A(cL 



A(c): 



Then, 



21 2 1 

Figure 2: Construction of the finite configuration 

v I I 

J~] c(i) = x Cp(i) + Cpp re fi Xi (») + Cpsuffi^ (« 



i=l 



i=l 



and 



^A(c)(i) = A(cp){i) + ^A(c pPrejix[ )(i) + ^Z yl ( 2 pSuffix i )(2 



i=l 



i=l 



z/3 <2«s 

We obtain the following inequalities: 

< < ra + 2/s and 



<2«s 



iG2 



Suppose that a < j3. If we choose an integer x such that x > -Jz^, we 
are in a contradiction with the previous inequalities. 

Similarly, if a > (3 then the contradiction is obtained by choosing x > 



Als 
a-0- 

Hence, if A is not PNC, it is not FNC either. 



□ 
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Proposition 2 A cellular automaton A is NC if and only if it is FNC. 



Proof. We first show that if A is not FNC then it is not NC either. Let A 
be a cellular automaton with a neighborhood of radius I. If A is not FNC 
then there exist an integer 7 and a configuration c in Cp such that 

7 7+i 

5>«=I>(0^ £ A(cm = j2Mcm. 

jgZ i=— 7 i=— -7— Z igZ 

Let {i^n} be a sequence of windows of size tending to the infinity, centered 
on the non-quiescent part of c. As the non-quiescent part of c is finite, there 
exists k such that, for all n > k, 

7 7+2 

Mn(c) = C(») A{c){i) = fl n (A(c)). 

i=— 7 i=— 7— 2 

Mn( c ) 

Hence, lim — 7- rrTT ^ 1, so A is not NC. 

Now assume that A is not NC. Then there exists a configuration c such 
that — does not converge to 1. Now let us call frontier of the window 

Vn(A{c)) 

the stripe of width the diameter of the neighborhood at the window border. 
Let us now count £ n (c) the number of non-zero cells of c in this frontier. 

First case: Va, 3n, pL n < an. Thus there exist consecutive subsequences 
of unbounded size in which fj, n (c) is constant. Let us consider such a sub- 
sequence of length greater than the neighborhood. The frontier of the last 
window is thus full of zeros. Hence we can "cut off" the internal part, sur- 
round it by O's and thus obtain a finite configuration. As A is not NC, it is 
not FNC either. 

Second case: 3a, Vn, fi n > an. The idea is to "cut off" the windows to 
form periodic configurations. We can bound the perturbation on ^i n (A(x) by 
the number of non-zero cells in the frontier multiplied by the maximal state (a 
constant). We then remark that with the hypothesis above, £ n (c) = o(fi n (c)). 
Thus for n sufficiently big, the unbalance of the ratio ^^ c ^ is transmitted 
to the periodic configuration obtained. Therefore A is not PNC and so not 
FNC. □ 



3 A necessary and sufficient condition 

We present below a necessary and sufficient condition for number-conservation 
Our condition is rather complicated to construct. For the sake of simplicity, 
we first present it in dimension 2 with a square neighborhood 2x2. Then 
we give the condition for neighborhood n x m (proof sketched). 
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Proposition 3 Let A = (2, N, Q, f) be a 2D cellular automaton with N 



(0,0) 

(o,-i) 

only if 



(0,1) 
(1,-1) 



and thus f : Q 4 — »■ Q, A is number-conserving if and 





c 



b 



c d 




a 





a b 



Proof. We first show that the condition is necessary. Assume that A is NC, 
let us denote by \c the sum of the states in a finite configuration c of A. As 
A is in particular FNC (see Prop. Ill) 



/ 



\A 





a 6 

c d 





a + b + c + d. 



Moreover, 



a + b 



d = f 



+.f 
+f 




a 



a b 

a b 
c d 



+ f 



+ f 




b 

c d 




+ / 



+ / 



c 


a 



+ / 



d 


b 
d 



(1) 



Replacing a by yields 

b+c+d = f 

+f 




b 



+ / 



d 




+ / 




b 

c d 




+ / 



+ / 



c 


b 
d 



+ .f 



+ .f 




c 

b 

c d 



(2) 



Subtracting 
/ 



from (H) we obtain: 

= a +f 
-f 




b 




a 





c 



a b 



+ f 
~f 



b 

c d 

a 

c 



(3) 
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Setting b = in (g), we have also: 



d = f 



-f 




c 



C G? 



/ 



/ 



c 


c d 




+ / 




d 



Subtracting (|j) from (§) we obtain: 



b 

c d 



+f 
-f 



We replace this term in (||) 



/ 



a b 
c d 



= a + b 
-I 

If A is FNC then 



-/ 


a 





d 



b 



+ / 

-/ 





c d 

b 

d 




d 





c d 



tf 




-/ 


b 




b 





a b 



(4) 



(5) 



a 
c 





c 



(6) 



V 





b 

d 





b + d. 



After similar calculi as above with the substitution of b by 0, we obtain 



b 
d 



= b + f 

-I 




d 

b 
d 




d 



-/ 




b 



-f 




b 



(7) 



By replacing this term in (||), we obtain exactly the formula of our Propo- 
sition. 

Let us prove now that this condition implies PNC. Then by Prop. [|, 
the theorem is proved. Consider a periodic configuration. Now sum this 
condition on the period; the terms [/(. ..)—/(...)] simplify and we get 
exactly the condition for PNC hence for NC. □ 

We now give the condition for more general neighborhood in dimension 2. 
We assert that an analogous formula can be obtained for higher dimension. 
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Theorem 1 Let A = (2, N, Q, f) be a 2D cellular automaton with 



N 



and thus f : Q r 



( (0,0) (0,1) ... (0,n-l) \ 

(1,0) (1,1) ... (l,n-l) 

V(m-1,0) (m-1,1) ... (m-l,n-l) J 
— ► Q. ^4 is number- conserving if and only if 



( zi,i a;i,2 

^2,1 ^2,2 



^2,n 



^1,1 



\ %m,l •Em.,2 ■ ■ ■ %m,n / 



I ... 



+ EE/ 

i=0 j=0 







V 



£1,2 
£2,2 



£l,n— j — 1 
£2,n-j-l 

%m—i,n—j — l J 



+ EE/ 

i=0 j=0 



/ ... 1 
... 

S v 

j 




>i+l 



#2,1 
£3,1 



£m— i— 1,1 



2,n— j 
£3,71— j 



%m—i—l,n—j J 



( o 



m— 1 n— 1 

E E / 

*=0 j=0(»+j>0) 



1 



£1,1 

£2,1 



x 2,n-j 
^m—i,n—j 
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'+1 

£2,2 
^3,2 

•^m—i— 1,2 



^2 5 7i— j-i 

•^3,n— j— 1 
^m— i — l,n— j — 1 / 



Proof. We give only the sketch of the proof. In order to show that the 
condition is necessary, we start from the configuration 





















£1,1 


£1,2 • 


2*1, n 








£2,1 


£2,2 ■ 






















•^771,1 


£m.2 ■ 























As A is FNC, we have £i,i+£i,2 + - • ■+x ni!n = sum of states obtained in each 
cell after the application of A, that is the sum of all /(...); we substitute 
in this equation x\^\ by 0, then x\$ by 0, and so on, up to xi^ n . Then we 
obtain an analogous formula to (||) 



/ ^1,1 

£2,1 



£1,2 
£2,2 



Xl, 
£2,' 



( £1,2 £1,3 
£2,2 £2,3 



£l,l + £1.2 + • • ■ + £l,n 



£l,-, 
£2,-, 







+ 



\ £m,2 £l,2 ■ ■ ■ X„ hn / 

Then, in order to obtain the final formula, starting from the configuration 





















£1,2 


£1,3 • 










£2,2 


£2,3 ■ 






















£m,2 


£m,3 
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we establish the equation for 



/ x l,2 
X2,2 



%l,3 
X2,3 



X2,: 



\ 





\ X m ^2 ^m,3 • ■ ■ *£m,n / 

Then, we substitute by 0, then x\^, and so on, up to xi iH . 

The condition is sufficient: for periodic configurations, when calculating 
the sum of states on the period, every /-term is present two times. If it is 
present in the sum of /-terms with X\^\ (resp. #2,2) as first non-zero element, 
then it is present in the sum of /-terms with X\^ (resp. £2,1) as first non-zero 
element. We so have as sum of states on a period Yli^i Sj=i x i,j- Hence 
this condition implies PNC and by Prop, [I] the theorem is proved. □ 



Theorem 2 The property "number-conserving" is decidable in ©(slogs) for 
cellular automata of any dimension, where s denotes the size of the transition 
table. 

Proof. Remark that the number of terms of our condition in Th. [l] is propor- 
tional to the number of neighbors. Thus the set of conditions to be checked 
is of size proportional to the size of the local transition rule of the cellular 
automata, which is upper bounded by the size of the cellular automaton in 
the chosen representation. Access to the table costs at most O(logs). □ 

Remark that if cellular automata are given by a program computing their 
rule and not by the rule exhaustively, then our theorem does not hold. This 
drawback is often present in cellular automata theory. But if we consider 
that access to the transition table is performed in constant time, then the 
algorithm is linear. 



4 Dynamics 

NCA are often used for the simulation of real phenomena such as highway 
traffic flow, fluid dynamics and so on. Therefore, once we have established 
that our model is a NCA (by using the decision procedures given in the first 
part of the paper, for example), it is very important to forecast as much 
as possible its dynamics and, whenever possible, correlate experimental evi- 
dence (such as pattern growth) with other interesting theoretical properties. 
This is the matter of this second part of the paper. 



4.1 Definitions and background 

A dynamical system (X, F) consists of a compact metric space X and a con- 
tinuous self-map F. Denote d the metric on X and F n the n-fold composition 
of F with itself. 
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We first recall some classical definitions. A point x G X is an equiconti- 
nuity point for F if for any e > there exists 5 > such that for any y G X, 
d(x,y) < 5 implies that Vt G N, d(F t (x) , F t (y)) < e. If X is perfect and 
all of its points is an equicontinuity point for F then F is equicontinuous. 
A dynamical system (X, F) is sensitive to initial conditions if there exists 
e > such that for any x G X and any 6 > 0, < d(x, y) < 5 implies that 
there exists t G N such that d(F t (x) , F 1 (y)) > e. We recall that in a perfect 
space any system with no equicontinuity points is sensitive to initial condi- 
tions and vice- versa. (X, F) is expansive if there exists e > such that for 
any x,y G X, d(x,y) > implies that 3t G N such that d(F t (x) , F t (y)) > e. 
If for any pair of non-empty open set U, V there exists n G N such that 
F n (U) n V 7^ then F is transitive. A point x is periodic if there exists 
n G N such that F n (x) = x. F is regular if periodic points are dense in X. 

Properties like sensitivity, expansivity, transitivity and regularity are of- 
ten used as basic components of chaotic behavior. For example, Devaney, in 
a popular book |J, defines as chaotic all systems that are transitive, regular 
and sensitive to initial conditions. 

We are interested in the study of dynamical systems in symbolic spaces. 
S z is endowed with the product topology (also called Cantor topology) of a 
countable product of discrete space. For the sake of simplicity, we will study 
cellular automata in dimension k = 1, the generalization of results to higher 
dimensions is straightforward for most of the results; Propositions § to [H] 
require deeper adaptations although no new idea is needed. 

The metric d on S* defined as 

Vx,yeS Z , d(x,y)=2- n 

where n = inf {i > 0,x(i) ^ y(i) or x(—i) ^ y(—i)} induces the product 
topology on S z . Denote by S* the set of finite words on S. For w G S* , 
I (w) denotes the length of w. For t G Z,w G S*, the sets 

[w] t = [x G S z , x(t) = w(0), ... ,x(t + £ H) = w(£ (w))\ 

are called cylinders and form a basis for the product topology on 5 Z . The 
shift map a is often used as a paradigmatic example of chaotic symbolic 
system. It is defined as Vc G S z , Vi G Z, o~(c)(i) = c(i + 1). Cellular 
automata are exactly those continuous maps from to S z that commute 



with the shift (i.e., F o a = a o F) \14\. 



4.2 Classification of dynamical behavior 

In spite of their very simple definition cellular automata may have complex 
dynamical behavior (see for example |^8|). The classification of these dy- 
namical behaviors is known as the classification problem is one of the oldest, 
hardest and for this reason fascinating open problem in cellular automata 
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theory. During the years many authors have proposed partial solutions (be- 
cause of their undecidability) . Each classification focuses on a particular 
aspect of cellular automata theory. For example, if one is interested in sim- 
ulations of real life phenomena like car traffic (e.g. using number- conserving 
cellular automata that we call NCA for short), the study of dynamics on fi- 
nite patterns deserves a particular importance. In this context a well-known 
classification of cellular automata with quiescent state is given in [||, In 
these papers cellular automata are classified according to their pattern di- 
vergence. Three classes of behavior are defined: 

Cy. ViG C F lim^oo i (F t (x)) = 

C2: Vi S Cf sup igN ^ (F*(x)) < 00 

C3: 3x G Cf sup t€N £ (F 4 (x)) = 00 

People familiar with the study of complex analytic transformations will im- 
mediately remark the strong analogies with Julia and Fatou sets of complex 
transformations |^]. 

This classification is especially pertinent for NCA since is necessarily a 
quiescent state. 

In the study of chaotic behavior one often finds that complex behavior 
is generated from "simple" local interactions. Again, using the metaphor 
of car traffic, traffic-jams and their consequences are generated by a "bad" 



interaction of a relatively small number of cars. In [18|, Kurka proposed 
a classification based on local behavior of cellular automata and increasing 
degrees of chaos. Kurka proposed the following classes: 

K\: equicontinuous cellular automata; 

K2: cellular automata with equicontinuity points but not equicontinuous; 
K3: cellular automata sensitive to initial conditions but not expansive; 
Ky. expansive cellular automata. 

The following result is adapted from Knudsen [17]. By using Proposition || 



one can express the properties defining (most of) Kurka classes in term of 
behavior on finite patterns. 

Proposition 4 F is transitive (resp. sensible to initial conditions, equi- 
continuous) w.r.t. finite configurations if and only if F is transitive (resp. 
sensible to initial conditions, equicontinuous) . 

The same result holds if in Proposition |4] we consider the set of spatial 
periodic configurations instead of finite configurations. In the rest of the 
section, we study the relations between dynamical behavior of NCA and the 
classifications of Kurka and Cattaneo. 
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Given x G Cf, let 



m(x) = min{i G Z, x(i) ^ 0} , M(x) = max{i € Z, ac(i) 7^ 0} and 

M(s) 

m(x) 

Remark that Vx G Cf, m(cr fc (2;)) = m(x) — k (the same holds if we put M 
in place of m). 

Proposition 5 T/ie are no expansive NCA. 

Proof. For dimension greater than 1 , this is a consequence of the well-known 
result of Shereshevsky |27y. Consider a one-dimensional NCA F. Suppose 
that F is expansive with expansive constant e. Let % G N be such that i < e. 
Construct a finite configuration Cj(j) = 1 if i = j, otherwise. Since F is 
expansive there exists t\ such that cZ(F fl (0), F' 1 ^)) > e. For the same 
reason there exists £2 such that d(F' 2 (0), F* 2 (c_j)) > e. Therefore there 
exists a time t > max{ii,i2} for which F*(cj) consists of at least two l's, 
violating FNC condition. □ 

Lemma 1 Let F be a NCA. For any x G Cf, there exists c x G N such that 
Vt G N, £(F*(x)) > c^. 

Proof. Let x £ Cf- Since each cell can contribute at most g — 1 to the sum 
s(x), it holds that Vt G N, I (F t (x)) > \^&\. □ 

The following proposition is a trivial consequence of Lemma [l[ 

Proposition 6 There are no NCA in class C\. 

Proposition 7 There are no NCA in class K\ fl C3. 

Proof. Let F be a cellular automaton in C3 and let c be such that 

lim t(F*{c)) =00 . (8) 

We have two cases: 

1. supm(F*(c)) = —00; 

teN 

2. sup m(F*(c)) = c G Z and sup M(F*(c)) = 00 . 

teN tefi 
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We consider only the former case, the proof for the latter is similar. For any 
5 > 0, let n G N be such that 4 < 5 < By (|8|), there exists a time t such 
that m(F*(c)) = -h > n. Therefore d(0, <r~ h (c)) < 5. Since F commutes 
with a, we have 

m{F\a- h {c))) = m{a- h {F\c))) = m{F l {c)) + h = 

which implies that 

d(F t ( ff -' l (0)),F t ( ( 7- ft ( C ))) = d(0,iiV- h (c))) = 1 • 

Since 5 is arbitrary, last inequality implies that is not an equicontinuity 
point. Since we are in a perfect space we conclude that F is not an equicon- 
tinuous system. □ 

A word W G £* is called blocking word if there exists an infinite sequence of 
words w n such that 

1. Vn G N, I (w n ) = 2h + l>r with h G N; 

2. Vc G (c_ fc:fc = W => Vt G iV \ {0} F*(W0_ i:l = ^) . 

In other words partitions the evolution diagram of F in two completely 
disconnected parts: perturbations made in one part are completely "blocked" 
by W. This notion is fundamental for characterizing cellular automata dy- 
namics in the presence of equicontinuity points as we can see in the next 
proposition. It will also be crucial in the proof of Proposition 11. 

Proposition 8 (|l[]) Any equicontinuity point has an occurrence of a block- 
ing word. Conversely if there exist blocking words, then any point with in- 
finitely many occurrences of a blocking word to the left and to the right (of 
0) is an equicontinuity point. 

Proposition 9 (fll]]) Let X be a a -transitive SFT. If F has an equicontinu- 
ity point then (F onto 4=> regular). 

The above proposition implies that surjective NCA in class K2 PI C3 are 
regular. 

Proposition 10 (0]) A cellular automaton F is equicontinuous and onto 
if and only if3p>0 such that Vc G S z , F p (c) = c. 

A consequence of the above proposition is that surjective NCA in class 
K\ (~l Ci are periodic and non-surjective ones are ultimately periodic. 

Proposition 11 There are no NCA in class Ki^Ci- 
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Proof. Since we are in a perfect space if a cellular automaton is not equiconti- 
nuous then there exists at least a configuration c which is not an equicontinu- 
ity point. By Proposition [| we can assume that c is finite. By the hypothesis, 
there exists a constant K such that Vi G N, £(F n (c)) < K. Assume that 
lim^oo m(F t (c)) = — oo (if this is not the case then lim^oo m(F t {c)) = +00 
and the rest of the proof will be perfectly similar). On the other hand, there 
exists at least an equicontinuity point and, by Proposition |§L there exists at 
least a blocking word W . Construct a configuration d such that 

(c(i) if m(c) <i< M(c) 
W{i) if 2M(c) < i < 2M(c) + 1 (b) 
otherwise . 

It is clear that lim I (F t (d)) = +00. □ 

For any cellular automaton F, denote Fp the restriction of F to Cp and 
Fp the restriction of F to Cp. 

Proposition 12 A NCA F is surjective if and only if its restriction to Fp 
is bijective. 



Proof. One implication is the well-known Moore-Myhill theorem (19 



For the converse, remark that if F is surjective every finite configuration has 
a pre-image. The FNC condition grants that this pre-image is finite. □ 

The following definition is essentially taken from @. 

Definition 4 A cellular automaton < S z , G > is said to be a generalized 
subshift on the set U C S z iff 

1. U is G positively invariant i.e., G(U) C U; 

2. there exists a mapping T, T' : S z -> N\{0} such thatMc £ U, G T{ - C \c) = 
a T '( c \c). 

The following lemma present peculiar properties of the shifting function 
of generalized subshifts on Cp. 

Lemma 2 Let F be a generalized subshift on Cf and let T,T' as in Defini- 
tion^. Then either Vx G Cf,T(x) = T'(x) or Vx G Cf,T(x) = —T'(x). 

Proof. By contraposition, suppose that there exist x, y G X such that 
T(x) = T'(x) and T(y) = -T'(y). Let n x = max <j< r(a .) t [F % {x)) and 
n y = max <j<T(?/) £ [F i (y)). Build the following finite configuration 

y(i + £ (y)) if - £ (y) - n y < i < -1 - n y 
Vi G Z, w(i) = ^ x(i — 1) if 1 + n x < i < £ (x) + n x 
otherwise . 
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Since lim^oo £ [F t (w)) = +00 there exists none Z\{0} such that F n (w) = 
a n (w) or F n (w) = (j- n {w). □ 

The following lemma links properties of the evolution of F with set properties 
of the maps Fp and Fp. 

Lemma 3 Let F be a transitive cellular automaton. Then Fp is bijective 
and Fp is surjective. 

Proof. It is well-known that transitivity implies surjectivity. By Moore- 
Myhill theorem, this last property implies that Fp is injective. By Proposi- 
tion [|, transitive cellular automata are transitive over Cp, This fact implies 
that Fp is surjective. Using Proposition [| again, we have that transitive 
cellular automata are transitive over Cp and hence Fp is surjective. □ 



Proposition 13 Generalized subshift cellular automata on Cp are transitive 
and regular. 

Proof. Let F be a generalized subshift cellular automaton on Cp and T, T' 
as in Definition [|. We are going to prove that F is transitive on Cp; by 
Proposition [| we have the first part of the thesis. For any cylinder set 
[w]t, let w the finite configuration such that Wi + t = Wi + t for < i < £{w) 
and elsewhere. Take two arbitrary cylinder sets [w]t and [z]i. Let = 
max <j< T (^) i (F l (w)) and n 2 = max <j< T (2) £ (F l (z)). Let m = T(w)T(z). 
Finally, choose k,q 6 Z such that q > n 5 + and t = km + q. Build a 
finite configuration c such that c = w © a krn+q (z), where © is the bitwise 
"or" operation. It is not difficult to see that F fcm+< ?(c) <E [z] and c G [wj; this 
completes the first part of the proof. 

Since Cp is dense in S*, for any configuration x G 5*, there exist a 
sequence {uj} igN C C^? such that lim^ootij = x. and let 

n Ui = max £ (F j (ui)) . 
0<j<T( Ui ) V ' 

For any «j, build a spatial periodic configuration 

y . = 2n n Xi 2n *i() 2n * iXi 2n *i . . . 2n ^Xi0 2n *> . . . 

such that (yi)k = (xi)k for m(iEj) < k < M{xi). Clearly yj is ultimately 
periodic for F. By Lemma || we have that it is in fact periodic. The fact 
that lim^oo yi = x closes the proof. □ 

From Proposition 13 and |12, 18], generalized subshift cellular automata on 
Cf are sensitive to initial conditions. 

When studying asymptotic behavior of a cellular automaton it is of fun- 
damental importance to understand the structure of the so called limit set 
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= n nG N-P n (5' z ). This set can be viewed as the set of configurations which 
have an infinite number of predecessors. It is well-known that is compact, 
non-empty and that it is the maximal (w.r.t. set inclusion) attracting set. 
In the sequel we will very interested in a special subset of the limit set: 
£If = n n€ fqF n (Cp) C Cp, that is the set of finite configurations which have 
an infinite number of predecessors (of finite type). One can easily see that 
for any cellular automaton this set is non-empty for it should contain at least 
one homogeneous configuration. In the particular case of NCA, the follow- 
ing easy proposition proves that the cardinality of ftp is always countably 
infinite. 

Proposition 14 For any NCA, O is an uncountable set while the cardinality 
ofVlF is countably infinite. 

Proof. Let c\ be the configuration that has a 1 in zero and elsewhere. 
Let h = \m(F^ci))\ (for a G Z, \a\ is the absolute value of a). Let = 
2h+i b0 2h+i for b e {0,1}. For any configuration x G {0, 1} Z , build the 
configuration c x = a Xf) a x , ...a x ... Note that c x G O. In particular if 
x G Cp (~1 {0, 1} then c x G tip. It is easy to see for x, y G {0, 1} , if x ^ y 
then c x ^ c y . □ 

A cellular automaton F is said to be a ultimately generalized subshift 
on U C S z is F is a generalized subshift on U and Vx G <S Z , there exists 
n = n x G N such that F Ux (x) G U. 

Proposition 15 Non-surjective NCA in K3 n C2 are ultimately generalized 
subshifts onQp. All surjective NCA in class K^nC2 are generalized subshifts 
on Cp. 

Proof. Consider a NCA F in K%C\C2, we prove that it is a generalized subshift 
on Uf C Cp. Since F is not equicontinuous for any finite configuration c 
either m(F t (c)) or M(F t (c)) is not a constant function. Since F G C2 there 
exists L G N such that Vi G N, £ (F t (c)) < L. Therefore there exist tx,t 2 < 
S L such that F tx {F t2 {c)) = a tl (F t2 (c)). Clearly, Vt F is positively invariant 
and, by the above remarks, F is a generalized subshift on it, for otherwise 
one can find a finite configuration in which m{ ) is constant. Moreover, since 
Vi G N, I {F\c)) < L, there exists t G N such that F\c) G Up. For the 
second part of the proof remark that surjective NCA are bijective on Cp (cf. 
Proposition |l^) then Of = Cp. □ 

Proposition 16 Let F be an ultimately generalized subshift on U C Cf with 
= co- Then F is sensitive to initial conditions. 
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Proof. We prove that they are sensitive to initial conditions on Cf, by 
Proposition ||] we have the thesis. 

Consider a finite configuration c and B$ the ball of radius 5 > centered 
on c. Choose an integer n such that ^ < 5. Let t G N be such that F'(c) G 
Of (by the hypothesis, t always exists). Let n c = max 0<i< y( c ) +i I (F*(c)) . 
Choose £ G fi_p. By Lemma without loss of generality we can assume that 
the patterns generated by c and x move to the left. Now, choose h > t such 
that m(F h (c))+n c < 0. Let n x = max < i<T ( 3; ) I (F l (c)) . Choose k such that 
kT{x) > h + n x + n and m(F kT( - x \x)) = 0. Finally, build the configuration 
w = c © <7 fcT ( x )(a;), where © is the usual bit-wise "or" operation. Is is not 
difficult to see that d(c,w) < 5 and d{F kT{ ~ x \c), F kT W(w)) > 1. Since c and 
8 were chosen arbitrarily we conclude that F is sensible to initial conditions. 
□ 

The following corollary is an immediate consequence of Propositions [u] 
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Corollary 1 NCA of class C2 that are ultimately subshift on U C Cf are 
sensitive to initial conditions. 

4.3 Examples 

This section proposes several examples which illustrate representative be- 
haviors in each class. 

Example 1 (Class K\ n C2 is not empty) 

K\ n C2 contains the identity cellular automaton.d 



Example 2 (Class K3 n C2 is not empty) 

Consider the NCA with the following local rule f on the alphabet S = {0, 1}: 



(x,y,z) 


000 


001 


010 


011 


100 


101 


110 


111 


f(x,y,z) 











1 


1 


1 





1 



Rule f is called "rule 184" in j^J and it often given an elementary example of 
car traffic modeling: symbols l's are cars moving on a highway (0's). When 
a 1 finds a on its right (that a free path on the highway) it moves one cell 
to the right. If the cell is not free then it waits. □ 



Example 3 (Class K2 (~l C3 is not empty) 

Consider the following local rule f on the alphabet S = {0, 1, 2}: 



(*) 


2a 


lb 


12 


0a 


m 


2 


b 


1 
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where a G A and b G {0, 1}. Let F be the global rule induced by f . It is clear 
that W = 2 is a blocking word. Define three configurations c, d, e respectively 
as c(0) = 2, c(l) = 1, otherwise; d(0) = 1, d(l) = 2, otherwise; and 

e(— 1) = 1, c(0) = 0, c(l) = 2, otherwise. Since sup£(F t (c)) = oo and 

ten 

F^ 1 (d) = {d, e}, by Moore-Myhill theorem JI^ ; 2i], F is not surjective and 
hence not regular. Finally, since is not an equicontinuity point we have 
thatF G K 2 nC 3 - □ 



Cellular automata in class K 3 n C3 seems to have a more rich variety 
of behaviors than other classes; as a consequence is much more difficult to 
have clear-cut situations: we can have system of particles with delayed prop- 
agation (Example ||) and others in which particles cross each other without 
interaction (Example ||) . We also underline that it is exactly in this class that 
reversible NCA capable of universal computation have been found 21]. 



Example 4 (Class K3 n C3 is not empty) 

Consider the NCA with the following local rule f on the alphabet S 
{0,1,2}: 



ix,y,z) 


00b 


002 


01a 


02a 


10a 


11a 


f(x,y,z) 





2 


a 





1 


a 



and 



{x,y,z) 


12a 


206 


202 


21a 


22a 


f{x,y,z) 


1 





2 


a 


2 



One can verify that for any finite configuration c there exists a time i such 
that for any t > t, F l (c) is such that all symbols 1 (if any) are on the 
positive coordinates, while symbols of type 2 are on the negative part. 1 's 
and 2 's propagate in opposite directions. Normally type 2 symbols propagate 
at speed 1 to the left (one unit displacement per time unit). This speed 
reduces of one half if 2 encounters a 1 like in 102; in a sense type 1 's symbols 
have precedence over 2 's. Therefore for any e > choose n G N such that 
= < e < • Construct the following configuration c as follows: 



c(i) ifm(c) < i < M(c) 
Mi G Z, c(i) = { 2 i = max{n, t} 
otherwise . 

Then there exists h < 2max{ra, t} such that d(F h (c) , F h (c)) 
= 1. From the fact that d(c,c) < e and Proposition we have that F is 
sensitive to initial conditions. Note that F is not surjective since its restric- 
tion to configurations on the alphabet {0, 1} coincides with elementary rule 
184 of Example ^ which is not surjective. □ 
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Example 5 Consider a system which describes the interaction of two par- 
ticles (denoted by 1 and 2) in a neutral media (symbols 0). Particles of type 
1 move left to right while type 2 move in the opposite direction. A symbol 3 
denotes the presence of a particle of type 2 and one of type 1 on the same 
site. These elementary rules are formalized by the following local rule: 





Oxa 


lxa 


2xa 


3xa 


f{x,y,z) 





1 





1 


(x,y,z) 


Oxb 


lxb 


2xb 


3xb 


fix,y,z) 


2 


3 


2 


3 



where a G {0, 1}, b € {2, 3} and x <E {1, 2, 3, 4} . 

It is easy to see that F is injective and hence its restriction to spatial 
periodic configuration is injective (see The last property implies that 

F is regular. F is also transitive. In fact, for any finite configuration x = 
X-k ■ ■ ■ %-k+h (k,h € N) construct a configuration y such that for —k — h < 
i < —k, y(i) = 1 or 3 if x(i + h) = 1, y(i) = otherwise. Moreover 
for —k + h < i < —k + 2h let y(i) = 2 if x(i — h) = 2 or 3. Remark 
that F h (y) = x. For any pair of cylinders C\ = C(w_fe, . . . ,W-k+h) an d 
C2 = C(z-. g , . . . , Z- g+n ) (k,h,g,n £ N). Choose c£C\. Let t be such that 
m(F t (c)) < -g + h and M(F t (c)) > -g + n + h. Define d such that for 
—g < i < —g + n, d(i) = z%. Moreover, for m(F t (c)) < i < m(F t (c)) + h 
let d(i) = 1 if Wi_ m{F t {c)) = 1 or 3; for M{F t {c)) -h <i < M(F*(c)) let 
d{i) = 2 if Wi_M(F t {c)) = 2 or 3. Finally, d{i) = in all other cases. By 
construction d G C2. Using previous remark one can see that F t {d)(j) = Wj 
for —k<j<—k + h, that is to say that F t (d) € C\. □ 

Table |] summarizes the situation. 





Ci 


c 2 


c 3 




(Prop. |6|) 


Example |l 




(Prop. 0) 


K 2 


(Prop. |6|) 


(Prop. 11 


) 


Example |3] 


K 3 


(Prop. |6|) 


Example || 


Example £| 


K A 


(Prop. |) 


(Prop. |) 


(Prop. |) 



Table 1: NCA in Kurka's and Cattaneo's classifications. 



5 Conclusions and open problems 

The paper discusses three formulations of the notion of being number- conser- 
ving which have arisen in recent works on the study of physical phenomena 
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ruled by conservation laws of mass or energy. We have proved that all these 
formulations are equivalent and we have given another equivalent definition 
in terms of density. 

Moreover, we have given a linear time algorithm for deciding if a CA 
(of any dimension and any state set) is number conserving. We stress that 
the linearity of the algorithm is something rather surprising since most of 
the interesting structural or dynamical properties of CA are known to be 
undecidable. 

The strong constraints imposed by number-conservation simplifies the 
study of asymptotic evolutions. Results can be conveniently interpreted in 
terms of particles system. Class K\ n C3 is composed by systems with "non- 
moving" particles. In class K3 D C2 we find systems in which particles move 
along only one direction with identical speed (of course speed may vary 
from system to system). Systems in class Ki n C3 allow particles in many 
different directions (possibly opposite) with the peculiarity that particles 
have blocking collisions. Finally, systems in class K3 n C3 have the same 
properties as those in class if 2 H C3 but with non-blocking collisions. 

Although many aspects have been understood there are still many inter- 
esting open problems that we think they worth further study. For example, 
we have already said that surjectivity property is undecidable for CA in 
general. But maybe it turns out to be decidable for the sub-class of NCA. 
Surjectivity plays a prominent role in the study of CA dynamics since it is a 
necessary condition (for example) for many popular definitions of chaotic dy- 
namics. In Section ^2] we have found that, with exception for class K3 n C3, 
surjectivity is equivalent to the property of having a dense set of periodic 
points (again, this last property is one of the component of Devaney's defi- 
nition of chaos). For class K3 n C3 we do not know if this property holds or 
not. 
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